A parametric approach to list decoding of 
Reed-Solomon codes using interpolation 

Mortuza Ali and Margreta Kuijper* 
April 15, 2011 

Abstract 

In this paper we present a minimal list decoding algorithm for Reed-Solomon (RS) codes. 
Minimal list decoding for a code C refers to list decoding with radius L, where L is the minimum 
of the distances between the received word r and any codeword in C. We consider the problem 
of determining the value of L as well as determining all the codewords at distance L. Our 
approach involves a parametrization of interpolating polynomials of a minimal Grobner basis 
G. We present two efficient ways to compute G. We also show that so-called re-encoding can 
be used to further reduce the complexity. We then demonstrate how our parametric approach 
can be solved by a computationally feasible rational curve fitting solution from a recent paper 
by Wu. Besides, we present an algorithm to compute the minimum multiplicity as well as the 
optimal values of the parameters associated with this multiplicity which results in overall savings 
in both memory and computation. 

1 Introduction 

Reed-Solomon (RS) codes are important linear block codes that are of significant theoretical and 
practical interest. A (n, k) RS code C, defined over a finite field F, is a k dimensional subspace of 
the n dimensional space F™. For a message polynomial m(x) = too + m\x + • • • + mk~ix k ~ 1 , the 
encoding operation is to evaluate m(x) at x\,X2, ■ ■ ■ ,x n , where the x^s are n distinct elements of 
F. The rich algebraic properties and geometric structures of RS codes lead to the invention of a 
number of efficient decoding algorithms such as Sugiyama algorithm [55], Berlekamp-Massey (BM) 
algorithm [H[2U], and Welch-Berlekamp (WB) algorithm [33]. These classical decoding algorithms 
guarantee correct decoding as long as the number of errors is upper bounded by t = \_(d — 1)/2J , 
where d — n — k + 1 is the minimum distance of the code. 

In classical decoding, the error correcting radius of t = \_(d — 1)/2J originates from the requirement 
of unique decoding since for t > \_(d — 1)/2J multiple codewords within distance t from the received 
word r may exist. One way to circumvent this limitation is to increase the decoding radius beyond 
[(d — 1)/2J and allow the decoder to output a list of codewords rather than one single codeword. 
However, such list decoding is only feasible if there are few codewords in the list. In |9j Guruswami 
showed that for a code of relative distance S = d/n, any Hamming sphere of radius < n(l — \J\ — 5) 
around a received word r contains only a polynomial number of codewords. Therefore, a (n, k) RS 
code with d = n — k + 1 can be list decoded up to the error correcting radius of n — y/n(k — 1) which 
Guruswami named as the Johnson bound. 
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A list decoding algorithm was first discovered for low rate RS codes by Sudan [35] and later improved 
and extended for all rates by Guruswami and Sudan 10J. The Guruswami-Sudan algorithm can 
correct errors up to the Johnson bound n — \J n(k — 1). Given a received word r, the essential idea 
of the algorithm is to find all the polynomials m of degree less than k such that m(x i ) 7^ for at 
most t values of i G {1, 2, . . . , n}. The Guruswami-Sudan algorithm finds these polynomials in two 
steps: the interpolation step and the factorization step. In the interpolation step, it computes a 
bivariate polynomial Q(x, r) that passes through all the points (xi, ri), {X2, rs), . . . , (x n , r n ) with a 
prescribed multiplicity s satisfying a certain weighted degree constraint (see |10j for the definition 
of weighted degree). Then the bivariate polynomial Q(x,r) is factorized to find all the factors of the 
form r — m(x), where to is a polynomial of degree less than k. Now a polynomial m is a valid message 
polynomial if it is of degree less than k and m(xi) ^ for at most t values of i 6 {1,2,..., n}. 
The construction of Q(x,r) with the prescribed multiplicity and weighted degree constraint ensures 
that for all valid message polynomials to, r — m{x) appears as a factor of Q(x,r). Even though the 
algorithm may produce implausible polynomials, the total number of polynomials L in the list will 
satisfy the bound L < (s + 0.5)-^/n/(fc — 1), see [2"T] . 

The most computationally intensive operation in the Guruswami-Sudan algorithm is the construc- 
tion of the bivariate polynomial Q(x,r). Computation of Q(x, r) involves solving a system of 0(ns 2 ) 
homogeneous equations which using Gaussian elimination can be done in time cubic in the number 
of equations [27] . Clearly the algorithmic complexity of the interpolation step is dominated by the 
multiplicity s. Recently Wu [29j transformed the interpolation problem to a 'rational interpola- 
tion problem' which involves smaller multiplicity. Given the received word r, Wu's algorithm first 
computes the syndrome s of r followed by the computation of the error locator polynomial A and 
error correction polynomial B using the Berlekamp-Massey algorithm. Wu demonstrated that all 
valid error locator polynomials can be expressed as a parametrization of A and B. More specifi- 
cally, given a list decoding radius t, Wu's algorithm aims at finding all polynomials A and j3 such 
that A' = AA + f3B has at most t distinct roots. Wu showed that similar to the Guruswami-Sudan 
approach, this problem can be reduced to a curve fitting problem but with significantly smaller 
multiplicity. 

It may be observed that the set of all Q(x,r) £ ¥[x,r] passing through the points (cEj,r,), for 
i = 1,2, ... ,n, with multiplicity s is an ideal / s . From this observation several authors including 
Alckhnovich [2], Nielsen and tfoholdt [23 , Kuijper and Polderman [TB], O'Keeffe and Fitzpatrick [23] , 
and Lee and O'Sullivan [IS], formulated the interpolation step of the list decoding algorithm as the 
problem of finding the minimal weight polynomial from the ideal I s . Clearly the minimal weight 
polynomial will appear as the minimal polynomial in a minimal Grobner basis of I s computed with 
respect to the corresponding weighted term order. Lee and O'Sullivan also showed that the minimal 
polynomial in the ideal I s can be computed more efficiently from a minimal Grobner basis of a 
submodule of ¥[x] q for a sufficiently large qM Let ¥[x, r] q = {/ € ¥[x, r] |r-deg(/) < q}. Then 
¥[x, r] q can be viewed as a free module over ¥[x] q with a free basis 1, r, . . . , r 9_1 . Then the essential 
observation of Lee and O'Sullivan is that the minimal polynomial of I s can be constructed from 
the minimal Grobner basis of a submodule of ¥[x] q along with the free basis 1, r, . . . , r 9_1 , for large 
enough q. 

In this paper we employ the theory of minimal Grobner bases to perform minimal list decoding. 
Given the received word r, let L denote the value of djj(r, C) where 

d H (r, C) := min{d H (r, c)}. 

Our main objective is to determine the value of L as well as all codewords c which are at a distance L 
from the received word r. Clearly, if L is larger than the classical error correcting radius L(^ _ 1 )/2j , 

Here the integer q is not related to the size of the field. 
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the task is a list decoding operation. Our algorithm, unlike the Lee and O'Sullivan approach, starts 
with computing a minimal Grobner basis G of a submodulc of F[x] 2 , rather than ¥[x] q . We then 
demonstrate that all valid message polynomials can be extracted from a parametrization in terms of 
the elements of G. For computational feasibility, we show that this parametric approach, like Wu's 
algorithm, can be translated into a 'rational interpolation problem'. However, our approach has at 
least three features that distinguishes it from Wu's algorithm. Firstly, our parametric formulation 
of the problem of list decoding, without the detour of syndrome computation, is simpler than Wu's 
formulation. Secondly, while Wu's algorithm, for each valid A', resorts to Forney's formula to 
compute the error values, our algorithm immediately leads to a valid message polynomial. Finally, 
we provide an algorithm to compute the minimum multiplicity along with the optimal values of the 
associated parameters to be used in the rational interpolation step. Use of these optimal parameters 
in the rational interpolation step results in savings of both memory and computation as compared 
to Wu's algorithm. 

The organization of the rest of the paper is as follows. In Section [2j we briefly review the relevant 
theory of Grobner bases. In Section |3j we develop the theory and present the main algorithm 
along with two ways to compute the minimal Grobner basis. In this section we also explain how 
so-called re-encoding can be applied to the proposed approach. In Section [4j we translate the 
parametric approach into a 'rational interpolation problem' and present an efficient algorithm for 
the computation of the minimum value of the multiplicity and other parameters to be used in the 
rational interpolation step. We demonstrate that the use of these optimal values of the parameters 
results in less memory requirement as well as less computational requirement as compared to Wu's 
approach. Finally we conclude the paper in Section [5] 

2 Preliminaries 

The theory of Grobner bases for modules in ¥[x] q is generally recognized as a powerful conceptual 
and computational tool that plays a role similar to Euclidean division for modules in ¥[x]. More 
specifically, minimal Grobner bases prove themselves as an effective tool for various types of inter- 
polation problems. In recent papers jTHl HZ] this effectiveness was ascribed to a powerful property of 
minimal Grobner bases, explicitly identified as the 'Predictable Leading Monomial Property'. The 
proofs in this paper make use of this property. Before recalling the PLM property let us first recall 
some terminology on Grobner bases. 

Let ei, . . . , e q denote the unit vectors in ¥ q . The elements x a e; with i e {1, . . . , 1} and a € No are 
called monomials. Let n\,...,n q be nonnegative integers. In this paper we define the following 
two types of monomial orders: 

• The (n 1 , ■ ■ ■ , n q )- weighted term over position (top) order, defined as 

x a ei < ar Qj a + rij < (3 + nj or (a + i%i = (3 + rij and i < j). 

• The (m, • • • , n q )-weighted position over term (pot) order, defined as 

x a ei < x^ Bj :<^> i < j or (i = j and a + m < j3 + nj). 

Clearly, whatever order is chosen, every nonzero element / € ¥[x] q can be written uniquely as 

L 

i=i 

where L e N, the Cj's are nonzero elements of F for i — 1,...,L and the polynomial vectors 
Xi, . . . , Xl are monomials, ordered as X% > ■ ■ ■ > Xl- Using the terminology of [T] we define 
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• lm(/) := X\ as the leading monomial of / 

• lt(/) := c\X\ as the leading term of / 

• lc(/) := ci as the leading coefficient of / 

Writing X\ — x ai e^, where a% £ No and i\ £ {1, . . . , I}, we define 

• lpos(/) := i\ as the leading position of / 

• wdeg(/) := a.\ + as the weighted degree of /. 

Note that for zero weights n± = ■ ■ ■ = n q = the above orders coincide with the reflected versions 
of the standard TOP order and POT order, respectively, as introduced in the textbook pQ. 

Also note that, unlike with TOP, the introduction of weights does not change the POT ordering of 
monomials. In this paper, weighted POT order is needed only because we need the associated notion 
of 'weighted degree'. 

We now recall some basic definitions and results on Grobner bases, see PQ. Below we denote the 
submodule generated by a polynomial vector / by (/). 

Definition 2.1 Let F be a subset of¥[x] q . Then the submodule C{F), defined as 

C{F) := (lt(/) | / e F) 
is called the leading term submodule of F. 

Definition 2.2 Let M C ¥[x] q be a module and G C M. Then G is called a Grobner basis of M 

if 

C(G)=C{M). 

In order to define a concept of minimality we have the following definition. 

Definition 2.3 (pfl Def. 4.1.1]) Let ^ f € ¥[x] q and let F = {fx, . . . , f s } be a set of nonzero 
elements of¥[x] q . Let ctj 1 , . . . , Oj m £ Nq and f3j 1 , . . . , /3j be nonzero elements of¥, where 1 < m < s 
and 1 < ji < s for i = 1, . . . , m, such that 

1. lm(/) = x aj i lm(/j- 4 ) for i = 1, . . . ,m and 

2. lt(/) = fax** lt(f n ) + ■■■ + im a;«*» lt(/ Jm ). 
Define 

h:=f- (fi h x^f h +■■■+ P jm x a ^f jm ). 
Then we say that f reduces to h modulo F in one step and we write 

If f cannot be reduced modulo F, we say that f is minimal with respect to F. 

Lemma 2.4 ( 1, Lemma 4.1.3]) Let f, h and F be as in the above definition. If f — > h then h = 
or \m(h) < lm(/). 

Definition 2.5 (pQ) A Grobner basis G is called minimal if all its elements g are minimal with 
respect to G\{g}. 
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It is well known [TJ Exercise 4.1.9] that a minimal Grobner basis exists for any module in ¥[x] q and 
that all leading positions of its elements are different. In [THl [T7] another important property of a 
minimal Grobner basis is identified; the theorem below merely formulates a well known result. 

Theorem 2.6 ( |18| ) Let M be a submodule of¥[x] q with minimal Grobner basis G = {gi, . . . ,g m }- 
Then for any 7^ / G M , written as 

f = a\9i H V a m g m , (1) 

where ai, . . . , a m € ¥[x], we have 

M/) = max , (M a i) lm (5i))- ( 2 ) 

l<z<m;ai^tO 

The property outlined in the above theorem is called the Predictable Leading Monomial (PLM) 

property, as in [18]. Note that this property involves not only degree information (as in the 
'predictable degree property' first introduced in (5]) but also leading position information. Most 
importantly, the above theorem holds irrespective of which monomial orders top or pot is chosen, 
for a proof see [TH] . 

Clearly, in the above theorem m = rank (M) and all minimal Grobner bases of M must have 
rank (M) elements, no matter which monomial order is chosen. Furthermore, we have the following 
theorem. 

Theorem 2.7 Let ni,...,n q be nonnegative integers and let M be a module in ¥[x] q . Let G = 
{<?l ; ■ ■ ■ t 9m} be a minimal Grobner basis of M with respect to the (m, • • • ,n q )-weighted top order; 
denote := wdegg.; for i = 1,.. . ,m. Let G — {<?i, . . . ,g m } be a minimal Grobner basis of M with 
respect to the (rti, • • • , n q )- weighted pot order; denote ti := wdeg <?i for i = 1, . . . ,m. Then 

m m 

£a = X> (3) 

i=l i=l 



Proof We first prove the theorem for the case m = q. It follows easily from the fact that both 
G and G are bases for M (in a linear algebraic sense) that there exists a unimodular polynomial 
matrix U e¥[x] qxq such that 

col {#!, ...,g q } = [/col {51, . . .,g q }. 

Without restrictions we may assume that the leading positions within each Grobner basis are strictly 
increasing. Clearly it follows from the above equation that also 

V = UW, (4) 

where V = col {gi, . . . ,g g }diag {x ni , • • • , x n "} and W = col {§1, . . . , g g }diag {x ni , • • • , x 7lq } Since 
U is unimodular we must have deg det V = deg det W. Clearly deg det V — YliLi ^ an d deg det 
W = Y^iLi &i fr° m which ^ follows. Next, we prove the general case m < q. For this, we note that 
it follows immediately from Q that the maximum degree of all minors of V equals the maximum 
degree of all minors of W. On the other hand, the maximum degree of all minors of V clearly equals 
SfcLiA an d similarly the maximum degree of all minors of W equals X)I=i^i- The theorem now 
follows. □ 



We call the sum in (J3J) the (rij, • • • , n q )-weighted degree of M, denoted by wdeg (M). For zero 
weights rii = • • • = n q = the above result expresses that the sum of the degrees of a (reflected) 
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TOP minimal Grobner basis of a module M coincides with the sum of the degrees of a (reflected) 
POT minimal Grobner basis of M. This result is merely a reformulation of the well known fact that 
the McMillan degree of a row reduced polynomial matrix equals the sum of its row degrees, see [6]. 

Corollary 2.8 let M be a module in ¥[x] q . Let G — {g\, . . . ,g m } be a Grobner basis of M whose 
(ni, • • ■ , rig) -weighted top degrees add up to wdeg (M). Then G is a minimal Grobner basis of M 
with respect to the (tlx, ■ ■ ■ ,n q )-weighted top order. 

Proof Suppose that G is not minimal. Then there exists g G G that can be reduced modulo G\{g}. 
This implies that there exists a Grobner basis of M whose sum of weighted degrees is strictly less 
than wdeg (M), which contradicts the above theorem. □ 



3 Minimal list decoding through division 

Let us now consider a (n, k) RS code and a nonnegative integer t. The problem of 'list decoding up 
to t errors' is the following: 

List Decoding Problem: Given a received word (ri, • • • , r n ) G F n , find all polynomials m £ F[i] 
of degree < k such that 

m{xi) — ri for at least n — t values of i G {1, ... , n}. 
3.1 Main approach 

We introduce the following two polynomials in ¥[x]: 

n 

U(x) =J[(x - Xi ), (5) 
i=l 

and C as the Lagrange interpolating polynomial, i.e., the polynomial of least degree for which 

C(xi) = ri for all z G {1, . . . , n}. (6) 

Definition 3.1 Let r = (ri,-- - ,r n ) G F™. The interpolation module M(r) is given by the 
module in ¥[x] 2 that is spanned by the vectors gi := [II(x) 0] and g<i := [C(x) — 1]. 

Note that {51,52} is a minimal pot Grobner basis for M(r). The above defined interpolation module 
is crucial to our approach. With g 2 we associate the bivariate polynomial Q2(x, y) = C(x) — y; clearly 
Q2(xi, r^ — for all i G {1, . . . , n}. Similarly, with g\ we associate the polynomial Q\(x, y) = Tl(x); 
trivially Qi(xi, r^) = for all i G {1, ... ,71}. Now consider an arbitrary bivariate polynomial Q of 
the form Q(x,y) — N(x) — D(x)y for which Q(xi,ri) — for all i G {1, . . . ,n}. It can be shown, 
see [16], that [N — D] G M(r). Recall that list decoding up to t errors amounts to finding all 
polynomials m G ¥[x] of degree < k such that 

m(xi) — ri for all i G {1, . . . , n} except i = j±, . . . , ]l with L < t. 

In our context this amounts to looking for an interpolating bivariate polynomial Q of the form 
Q(x,y) = D(x)m(x) — D(x)y, where D(x) = \~\f = i{x — x^). Note that then indeed Q(xi,ri) = 
for all i G {1, . . . ,n}. Thus, to solve the above list decoding problem we are looking for particular 
vectors [N — D] G M(r) of weighted (0, k — l)-degree < t + k — 1, that satisfy 
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1 . N is a multiple of D and 

2. D has L distinct zeros in F, where L denotes deg D. 

In this paper we are interested in finding the smallest value L = (r, C) for which list decoding is 
possible as well as performing the associated list decoding. Thus we occupy ourselves with maximum 
likelihood list decoding. We have the following theorem. 

Theorem 3.2 Let r = (r l7 -- - ,r n ) G F" be a received word and let M(r) be the corresponding 
interpolation module. Let f = If*- 1 ' f^] € ¥[x] 2 be a vector in M(r) of weighted (0, k — \)-degree 
L that satisfies the following 3 requirements: 

1. lpos(/) = 2, 

2. is a multiple of f^ and 

3. there is no vector in M(r) of weighted (0, k — l)-degree < L that satisfies requirements 1) and 
2). 

Then 



is a message polynomial corresponding to a minimal error pattern of L — k + 1 errors. 



Proof From lpos(/) = 2 it follows immediately that deg m < k and deg f^ = L — k + 1. It 
remains to prove that f^ has L — k + 1 distinct zeros in F. Since / £ M(r) there exist polynomials 
a and j3 such that 

" n o 



f = [a P] 



C -1 



(7) 



Observe that a and f3 do not have a common factor, otherwise the weighted degree of / would not 
be minimal (requirement 3). From ^ it follows that all — f^C = is a multiple of f^ by 
requirement 2. As a result, all is a multiple of f^ 2 \ Since a and j3 = —f^ have no common factor 



it follows that II must be a multiple of f^ 2 ' 
the theorem. 



i.e., f^ has L — k + 1 distinct zeros in F, which proves 



□ 



Lemma 3.3 Let r = (n, ■ • • ,r n ) G F" be a received word and let M(r) be the corresponding in- 
terpolation module. Let {31,32} be a (0, k — l)-weighted top minimal Grobner basis for M(r) with 
lpos(<?2) = 2. Denote £1 :— wdegg! and £2 '■= w deg<?2- Let t be a nonnegative integer. Then a 
parametrization of all vectors f G ¥[x} 2 with lpos(/) = 2 and wdeg/ = t + k — 1 (with respect to the 
(0, k — l)-weighted top order) is given by 

f = ag x +bg 2 , 

where a G F[x] with deg a < t + k— 1 — £\ and b is a monic polynomial in ¥[x] of degree t + k — 1 — £2- 
In particular, there exist no such vectors f for t < £2 — k + 1 . 
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Algorithm 1 Minimal list decoding of (n, k) RS code 



Input: Received word r = (n, . . . , r n ) 

Output: A list of polynomials m of degree < k such that d^(c,r) is minimal, where c = 
(m(xi), . . .,m(x n )). 

1. Compute the polynomials II and L given by ^ and ([6| ; define the interpolation module 
M(r) := span {[n 0] , [C - 1]}. 

2. Compute a minimal Grobner basis G = {31, 32} of M(r) with respect to the (0, k — l)-weighted 
top monomial order, with lpos(g2) = 2. Denote l\ := wdeg 31 and £2 := wdeg 32; set j = 0. 

3 . Check requirement 2) of Theorem |3.2| for / = agi + bg2 , for all a € F [x] with deg a < £2 — £\ + j 
and for all monic b £ ¥[x] with deg b = j; write / = [f^ f^] ■ 

4. Whenever step 3) is successful, output all obtained quotient polynomials, i.e., polynomials m 
of the form m = — /^V/ • I R case s ^ e P 3) is not successful increase j by 1 and repeat step 3). 



Proof According to Theorem 2.6 {31,32} has the PLM property with respect to the (0, k — 1) 



weighted top order. The parametrization now follows immediately from this property. □ 

Together, the above lemma and theorem give rise to the heuristic list decoding Algorithm [l] 

An important feature of the above algorithm is that we use £2 = wdeg 32 to decide how many errors 
to decode. Indeed, it follows from the above lemma that it is not possible to perform list decoding 
for t < £2 — k + 1. We now present the main theorem of this section. 

Theorem 3.4 Let r = (ri,-- - ,r„) G F" be a received word and let M(r) be the corresponding 
interpolation module. Let {31, 32} be a (0, k — l)-weighted top minimal Grobner basis for M(r) with 

lpos(g2) = 2. Write 32 = g^ ■ Then Algorithm jij yields a list of all message polynomials m 

such that 

dn(c,r) is minimal, where c — (m(xx), . . . ,m(x n )). (8) 

In particular, in case there exists an error pattern with only < \_{n — k)/2\ errors, the list consists 
of only 

o (1) 

m = -%- (9) 

.92 



Proof Firstly, it follows immediately from Theorem |3.2| and Lemma |3.3| that any polynomial m 
that is output by Algorithm [l] has to have degree < k and satisfy Vice versa, if m is a polynomial 
of degree < k that satisfies Q then it follows from Lemma [3. 3| that it must be in the output list of 
Algorithm[T] Finally, let us assume that there are only < [_(n — k) /2J errors. This implies that there 
exists a vector / = [/M f^\ in M(r) with wdeg / < [(n - fc)/2j + k - 1 < (n + k - l)/2 that 



satisfies the requirements of Theorem 3.2 Because of Lemma [373] it follows that £2 < (n + k — l)/2. 
Now, since £\ + £2 = n + k — 1 by Theorem |2.7[ this implies that l\ > £2. As a result, a = in step 
3), so that step 4) immediately gives the unique solution for j = as ([9|. □ 

Our next example illustrates the classical decoding scenario, showing that Algorithm[l]is an extension 
of existing classical interpolation-based algorithms as in [THl [7] . 

Example 3.5 Consider the single-error correcting (7,5) RS code over GF(7). The message poly- 
nomial m(x) = 2x 2 + x + 3 is encoded as c = (m(0), m(l), . . . , m(6)) = (3, —1, —1, 3, —3, 2, —3). Let 
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the received word be r = (3, 2, —1, 3, —3, 2, —3). Thus an error occurred at locator position 1. The 
polynomials C and II are computed as C(x) — — 3x 6 — 3x 5 — 3x 4 — 3x 3 — x 2 — 2a; + 3 and II(x) = x 7 — x. 
Thus the module M(r) is spanned by the rows of the matrix 

x 7 — x 
-3x 6 - 3x 5 - 3a; 4 - 3a; 3 - x 2 - 2x + 3 -1 

A minimal Grobner basis {31,52} of M(t) with respect to the (0,4) -weighted top monomial order is 
computed as 



col {31,32} = 



-3x 6 - 3x 5 - 3a; 4 - 3a; 3 - x 2 - 2x + 3 -1 

-x + 1 



Thus, in the terminology of Theorem 3.4 we have 32 = 2a; 3 — x 2 + 2x — 3 and g^ = —x + 1. 
Applying Algorithm^^ we determine that is a multiple of and we recover 

m(x) — — '-J2) = 2x 2 + x + 3. 

92 

Let us now move on to an example of decoding beyond the classical error bound. Our approach is 
particularly feasible for the case that 6=1 and a is restricted to a constant, as illustrated in the 
next example. Note that the example is an instance of "one-step-ahead" list decoding [2"9"] . 



Example 3.6 Consider the single-error correcting (7,4) RS code over GF(7); let the message poly- 
nomial be m(x) = 2x 2 +x+3 which is encoded as c = (m(0), m(l), . . . , m(6)) = (3,-1,-1,3, —3, 2, —3) . 
Let the received word be r — (3, 2, —1, 3, 2, 2, —3) which differs from c at locations 1 and 4- The poly- 
nomials C and II are computed as C{x) = — x e — 2x 5 + x 4 — x 3 + 2x + 3 and H(x) — x 7 — x. The 
interpolation module M(r) is spanned by the rows of the matrix 



M(r) 



x 

~x 3 + 2x + 3 -1 



A minimal Grobner basis {31,32} of M(r) with respect to the (0, 3)-weighted top monomial ordering 
is computed as 



3 _ ^2 _l „ q _ 3x _ 1 
2x 4 + 2x A + x 2 - 3x + 2 x 2 + 2x - 3 



col {31,32}= o„..i , ... . , ..» ,2 



T/itts m this example l\ — £2 — 5, so that a is a constant. Applying Algorithm Q] we consider 
f = agi + 32 for a = 0, ... ,6. Writing f — [/W / 2 ] , we find that f^ divides /W J or a = 0, 2, 
and 4 7 giving a list of three message polynomials — we recover not only m(x) = 2x 2 + x + 3 (for 
a = 0), but also the message polynomials 3x 3 — 2.x 2 + 3x — 2 (for a = 2), and — 2x 3 — 2x 2 + 3x + 3 
(for a = A). 



3.2 Computation of gi and g 2 

There are various ways in which the required minimal Grobner basis {31,32} of the interpolation 
module M(r) can be computed. One obvious way is to simply run an existing computer algebra 
system such as Singular, specifying the required (0, k — l)-weighted top order. 

Because of the specific form of M(r) a more efficient way is to apply the Euclidean algorithm to the 
polynomials II and C. More specifically, we have the following algorithm. 
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Algorithm 2 Computation of gi and g 2 via Euclidean algorithm 



Input: Received word r = (7*1, . . . , r n ); polynomials II and C given by ([5) and 

Output: Polynomials <?i and <?2 in F[cc] 2 , such that {51,52} is a minimal Grobner basis of M(r) 

with respect to the (0, fc — l)-weighted top monomial order, with lpos(g2) = 2. 

1. Define polynomials ho, hi, to and ti in F[x] as 



ho 


to ' 




" n 


hi 






£ -1 



set j := 0. 

2. Check 

degt J+ i + k - 1 > deg h j+ i; (10) 

if NO, go to Step 3. If YES, define gi := [hj tj] and g 2 := [h j+x t j+ i] and STOP. 

3. Apply the Euclidean algorithm to hj and hj + i, yielding hj = qj+ihj+x+hj+2, where deg hj+2 < 
deg h j+ i. 

4. Write 



h j+ i 


tj+i 




' 


1 




hj tj 


hj +2 


tj+2 




1 


-Qj+i _ 




hj+i tj+i 



increase j by 1 and go back to Step 2. 



Theorem 3.7 Let r = (ri,-- - ,r„) G F™ be a received word and let M(r) be the corresponding 
interpolation module. Then Algorithm^ yields a (0, k — l)-weighted top minimal Grobner basis 
{51,52} for M(r) with lpos(.g 2 ) = 2. 



Proof Firstly we note that the matrix 



is unimodular, 
matrix 



1 

1 -qj+i 

i.e., has a polynomial inverse. It then follows that, at each step j, the rows of the 



■i j+ i t j+ i 

are a pot minimal Grobner basis for M(r) whose (0, k— l)-weightcd pot degrees add up to n + k— 1. 
By definition, with respect to the (0, k — l)-weighted top order both these row vectors have leading 
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position 1, until the stopping condition (10) is met. At this point the second row vector has leading 
position 2 and the sum of the (0, k — l)-weighted top degrees add up to n+k — 1. It now follows from 
Corollary 2.8 that the rows of the matrix (111 must be a (0, k — l)-weighted minimal top Grobner 
basis for M(r). □ 



Yet another alternative is to use an iterative method, interpolating the Xi's step by step for i = 
l,...,n. This method has the advantage that the Lagrange polynomial C does not need to be 
computed upfront. 

Theorem 3.8 Let r = (ri,-- - ,r„) G F™ be a received word and let M(r) be the corresponding 
interpolation module. Then Algorithm^ yields a (0, k — l)-weighted top minimal Grobner basis 
{51,52} for M(r) with \pos{g 2 ) = 2. 
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Algorithm 3 Computation of g\ and <?2 via iterative algorithm 



Input: Received word r = (7*1, . . . , r n ). 

Output: Polynomials g\ and <?2 in ^[a;] 2 , such that {<?i, 52} is a minimal Grobner basis of M(r) 
with respect to the (0, k — l)-weighted top monomial order, with lpos(g2) = 2. 



1. Initialize Lq 



1 and R := I e ¥ 



2x2. 



denote Ra 



Qj 



Nj -Dj 



€ ¥[x] 2x2 for 



j = 0, . . . , n. 

2. Process the received values rj iteratively for j = 1 to n as follows. For j = 1 to n do 

1. compute Tj := Qj—\(xj) — rjKj-i(xj) and Aj := Nj-i(xj) — TjDj-\[xj) 

2. define Rj := VjRj-i, where 



a, -r, 

X — Xj 



r J = o), 



• Vj := 







and Lj := Lj-i + 1 if Aj ^ and {Lj- X < (j + k - l)/2 



and Lj := Lj-i otherwise 



3. Define gi := [Q n - K n ] and g 2 := [N n 



Proof For j = 1, 



, n denote the interpolation module associated with r\ , 



by M (ri 



We show that the rows of Rj are a Grobner basis of M(y\, . . . , rj)of the required form for j = 1, . . . , n. 
We interpret Lj as the (0, k — l)-weighted top degree of the second row of Rj. Clearly this is true for 
j = l. Let us now proceed by induction and assume that this is true for j — 1 € {0, . . . , n — 1}. By 
definition of Vj and the induction assumption the rows of Rj are a Grobner basis for M(ri, . . . , r^). 
Also, by construction, their (0, k — l)-weighted top degrees add up to 1 more than the (0, k — 1)- 
weighted top degrees of Rj-i- Then, by induction, the (0, k — l)-weighted top degrees of Rj add up 
to j + k — 1 = wdeg (M(ri, . . . ,rj)). It then follows from Corollary 2.8 that the rows of Rj are a 
(0, k — l)-weighted top minimal Grobner basis for M(r l7 ...,Tj). Finally, by construction and the 
induction hypothesis, it is easily seen that the second row of Rj has leading position 2. This proves 
the theorem. □ 



3.3 The special case r = (y 1 , . . . , y n -k, 0, • • • ,0) 

In this subsection we pay special attention to the case that the received word r is of the form 
(yi, . . . , y n -k> 0, • • • , 0) € F". This comes about when so-called "re-encoding" is used in advance of 
RS decoding, see e.g., [TBI [12] . 

First we introduce the polynomial G € F[i] of degree k — 1 as 

n 

G:= J] (x-Xi). (12) 

i=n-fc+2 

Clearly, the polynomials II and C of the previous subsection can be written as 

n = Tl y G (13) 
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and 

C = C y G, (14) 
where 11^ and C y are in ¥[x]. The following lemma is straightforward. 

Lemma 3.9 Let (2/1, .. . , y n -k) G ¥ n ~ k , r = (j/i, . .. ,y„_ fc ,0, ■■■ ,0) 6 F" and Zet IT, £, G, IIj, and 
jC y oe defined as before. Let M(r) := span {[II 0] , [£ — 1]} as before and define M*(y) := 
span {[Ily 0] , [£ y — 1]}. T/ien the following two statements are equivalent: 

• {5i)ff2} *5 a minimal Grobner basis of M*(y) with respect to the unweighted top order, with 
lpos(# 2 ) = 2 

• {31,32} *s a minimal Grobner basis of M(t) with respect to the (0,k — \)-weighted top order, 
with lpos(</ 2 ) = 2, 



where gi 



(i) 



.,(2) 



and gi 



9 l 5; 



/ori= 1,2. 



Because of the above lemma it is now straightforward to modify Algorithm [T] into Algorithm |4j 
Algorithm 4 Minimal list decoding of (n, k) RS code for re-encoded received word 
Input: Received word y = (2/1, . . . , y n -k) in F n_fc . 

Output: A list of polynomials m of degree < k such that dn(c,r) is minimal, where c 
(m(xi), . . .,m{x n )) and r = (y u . . . ,y n -k,0 ,0) in F™. 



1. Compute the polynomials H y and C y given by (13) and (14) ; define the interpolation module 
M(y) := span{[n y 0],[C y - 1]}. 

2. Compute a minimal Grobner basis G = {ffi,<?2} of M(y) with respect to the unweighted top 
monomial order, with lpos((?2) = 2; set j = 0. 

3. Compute / = agi +652, for all a 6 F[x] with dega < £2 — 1\ +j and for all monic b E ¥[x] with 
deg6 = j; write / = [fW fW]. Check whether f^G is a multiple of f^ 2 \ where G is given 
by pi. 



4. Whenever step 3) is successful, output all obtained quotient polynomials, i.e., polynomials m 
of the form m = —f^G/ f^ 2 \ In case step 3) is not successful increase j by 1 and repeat step 3). 



Again the Euclidean algorithm can be used to compute g\ and 52; for this, Algorithm [2] should be 
initialized by 11^ and C y instead of II and C and the stopping criterion ( 10 1 should be replaced by 



deg t j+1 > deg hj+%, 



instead of ( 10 ) 



An alternative way to compute g\ and 172 is to employ an algorithm that processes the values of 
... , y n -k iteratively. For this, Algorithm[3]is modified into Algorithm[5]which essentially coincides 
with the well-known Welch-Bcrlekamp algorithm [28J , see also [TU [15] . 



4 Minimal list decoding through rational interpolation 

The most computationally intensive task in Algorithm [T] is Step 3. Recall that in Step 3, we need 
to determine all a and b of degree fci < £2 — @i + j an d k 2 — j such that /W is a multiple of f^ 2 \ A 
brute force approach may be to consider 



(2) 



(1) (2) , u (1) (2) 

a 9i 9i + b g\' g\ ! 
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Algorithm 5 Computation of gi and g 2 via iterative algorithm for re-encoded received word 



Input: Received word y = (y 1; . . . , y n -k) m F n_fe . 

Output: Polynomials g± and g 2 in F[x] 2 , such that {(71,(72} is a minimal Grobner basis of M(y) 
with respect to the unweighted top monomial order, with lpos(g2) = 2. 



1. Denote i? 7 



Qj -Kj 



for j = 0, . . . , n; initialize Lq := and 



i?o := 



cc - x„_ fe+ i 
1 



G F[af 



2. Process the received values j/j iteratively for j = 1 to n — k as follows. For j = 1 to n — k do 

1. compute Tj := Qj_i(xj) — rjKj-\(Xj) and Aj := Nj-i(xj) — TjDj-\(xj) 

2. define Rj := VjRj-i, where 



Aj -r, 

a; — a%- 



a; — 



and ij := Z^-i + 1 if Aj ^ and {Lj-\ < j/2 or Tj = 0), 
and Lj := otherwise 



3. Define gi := [<2„_fc - i^ n -fe] and g 2 := [iV n _ fc - -D n -fc]- 



and check for all polynomials a and & of bounded degree fci and fc 2 , respectively, whether /( 2 ) divides 
J^ 1 ). Clearly this approach is feasible only when both k\ and k 2 are small. For large values of k\ and 
k 2 , the computational complexity becomes prohibitively high, especially when the code is defined 
over a large field. Fortunately, Step 3 can be formulated as an algebraic curve fitting problem for 
which efficient polynomia l tim e algorithms exist. We explain this approach in the following. 
It follows from Theorem 3.2 that, in the context of Algorithm [l] is a multiple of if and 
only if /( 2 ) has t — £ 2 — k + 1 + j distinct roots. Therefore, an alternative approach to Step 3 is to 
determine all a and b of degree ki < t + k — l\ — 1 and k 2 — t + k — l 2 — 1, respectively, such that 

f ( - 2 \x)=a{x)gf\x) + b{x)gf\x) (15) 



(2) 

has t distinct roots. Now dividing both sides of (15) by g\ (x) we get 



5i 2) W Si 2) (z) 

Now let us define 



a(x)+b(xy-j^i. (16) 



9 2 \xt) 



for i = 1, • • • , n. 



Then Step 3 of Algorithm [T] can be formulated as the following rational interpolation problem. 

Rational Interpolation Problem: Given n points (xi, zi), (x 2 , z 2 ) ■ ■ ■ , (x n , z n ) and a non-negative 
integer t, determine all rational polynomials of the form z = a/b, with a and b of degree k\ and k 2l 
respectively, such that z passes through t of the n points (sci, Zi), {x 2l z 2 ) ■ ■ ■ , (x n , z n ). 
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This problem looks similar to the interpolation problem addressed by Guruswami and Sudan in jlOj . 
However, it is complicated by the fact that now we look for a rational solution rather than a polyno- 
mial solution. Recently, this rational interpolation problem has been addressed by Wu in [29]. For 
the sake of completeness we briefly describe Wu's formulation here. 



4.1 Wu's rational interpolation algorithm 



In line with the Guruswami-Sudan approach, Wu's algorithm first computes a bivariate polynomial 
Q(x, z), satisfying certain constraints, that passes through all the n points {x\, zi), (x2, z 2 ) • • • , (x n , z n ). 
Then the desired rational solutions z = a/b are obtained from the factorization of Q(x, z). Given 
the values of t, k\, and k%, let us determine the constraints that must be satisfied for the existence 
of such a Q(x, z). 

Let us define the (1, w) weighted degree of a bivariate polynomial Q(x, z) = ^2 



(i,j)£l u hJ 



wdeg x w Q(x, z) — max {i + jw}. 



(17) 



Let w := fci — k 2l p :— wdeg^ l w ^Q(x,z), and M := wdeg ± Q(x, z). Clearly deg 01 Q(x, z) is the 
z-degree of Q{x, z). Now if z = a/b passes through t points with multiplicity s then the polynomial 
b(x) M Q(x,a(x)/b(x)) must have ts roots. On the other hand, b(x)z — a(x) will be a factor of Q(x, z) 
if b(x) M Q(x,a(x)/b(x)) is identically zero. In turn, b(x) M Q(x, a(x)/b(x)) will be identically zero 
if it has more roots than its degree. Now the degree of b(x) M Q(x, a(x)/b(x)) is at most p + Mki 
Therefore, a necessary condition that must be satisfied is 



p + Mk 2 < ts. 



(18) 



On the other hand, a necessary condition for the existence of Q(x, z) passing through the n points 
with multiplicity s is that its (u, v)-th Hasse derivatives at all the n points are zero for all u + v < s. 
Thus the requirement that (xi,Zi) be a zero of Q(x,z) with multiplicity s, for all i = 1,2, ... ,n, 
leads to N constraints in the form of N homogeneous equations where 



N = ns(s + l)/2 



(19) 



and unknown variables are the coefficients of Q(x, z). A nonzero solution to the system of homoge- 
neous equations is guaranteed to exist if the number of equations is less than the number of unknowns. 
Now the number of coefficients in Q(x, z) with wdegj w Q{x, z) = p and wdeg 1 Q{x, z) = M is 



w 



U=(p+ 1)(M + 1) - —M(M + 1) 



(20) 



Therefore, a sufficient condition for the existence of a Q(x, z), passing through all the n points with 
multiplicity s, is 



(P + 1P + 1)-^(M + 1)>^. 



(21) 



Wu, in [29], has proposed suitable choices for the values of s, M, and p satisfying (18) and (21) as 

t(n-k + l-t) 



t 2 - n(2t -(n-k+l)) 



M = 
P = 



st 



2t-(n-k+ 1) 
ts - Mk 2 - 1. 



(22) 

(23) 
(24) 
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For more details on Wu's algorithm see 



It is worth noting that the multiplicity s, computed 



using (22), is not minimal. Although Wu suggested to first compute s according to (22) and then 



greedily minimize it subject to a certain constraint, he did not give any explicit algorithm to compute 
the minimal value of s. More importantly, given the minimum s, the values of M and p computed 
in (23) and (24) are not necessarily optimal. In the next section, we present an algorithm that 



computes the minimum value of s as well as the associated optimal values of M and p. 



4.2 Optimizing the integer parameters 

Given feasible values of s, M, and p, the rational interpolation step involves two steps: (1) con- 
struction of Q(x, z) and (2) factorization of Q(x, z). The best known algorithm for the construction 
of the interpolating polynomial Q(x,z) is the Kotter algorithm [11]. The Kotter algorithm has a 
complexity of 0(MN 2 ) [TJ], where N is the number of constraints as defined in ( 19 ). More precisely, 



it has memory complexity of 0{MU) and time complexity of O(NMU) [8:, where U is the number 
of coefficients in Q(x,z) as defined in (20) and M is the z-degree of the interpolating polynomial 



Q(x, z). On the other hand, the rational factorization step can be done in time 0(n 3 / 2 s 7 / 2 ) using 



Wu's rational factorization procedure [29]. As analyzed in sub-section 4.4 it is the Kotter algorithm 



that dominates the overall memory and computational complexity of the proposed, as well as Wu's, 
list decoding procedure. Therefore, to reduce the complexity of the Kotter algorithm, we take the 
following two step strategy. In the first step, we derive an explicit method to determine the minimum 
value of s for which there exist some M and p satisfying (18) and (21). Once the minimum multi- 



plicity is determined, N becomes fixed. Then in the second step, we compute the optimal values of 
M and p such that MU is minimized. 

The constraint ( 18 1 can be geometrically interpreted as follows. Assume that t and s are fixed. With 



the requirement that all the values involved in (18) are non-negative integers, all feasible values of p 



and M must be on or below the line L defined by the equation 

p + Mk 2 = ts - 1. 



(25) 



On the other hand, the constraint (21) requires that all feasible values of p and M are above the 
curve C defined by the equation 



(, + l)(M + l)^M(M + l)^. 



(26) 



Therefore, a necessary condition for the existence of a feasible solution satisfying both the con- 



straints (18) and (21) is that L intersects C at two different points (Mi,pi) and (M%, P2) on the real 



plane. Now solving (25) and (26) for M we get 

(ts - jfe ) ± V (ts - fco) 2 - 4(N - ts)k 



M = 



2k 



(27) 



where k = (k\ + k 2 )/2. According to Algorithm [T] while correcting t — l\ — k + 1 + j errors, we 
have k\ — £2 ~ £1 + j and k 2 — j- Using l\ + £2 = n + k — 1, we get ko = (t — to) where to = d/2. 
Substituting fc = (t — to) in (27) we get 



M 



(ts-t + t ) ± y/(ts - t + t ) 2 - A(N - ts)(t - t ) 
2(t - to) ' 



(28) 



It follows from ( 28 1 that the value of M and thus the choice of s is independent of k\ and &2 ■ Now 



for a fixed s, it can easily be verified if L and C intersect at two different points on the real plane 
by checking whether 

(ts~t + tof > A(N -ts)(t-t ). (29) 
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— 
— 



ns(s + 1) 




(p + l)(M + l)-|M (Af+ l) = 2(^) 




z-degree M 
00 

Figure 1: Consider correcting t = 7 errors in the decoding of a (15,5) RS code over GF(16) when 
fci = 2 and &2 = 1; (a) With s = 1, as the line L does not intersect the curve C, no feasible values for 
M and p exist; (b) With s = 7, the line L intersects the curve C at two different points (Mi = 14, pi) 
and (Ma = 17.67, ^2)- Thus M* — 15, 16, 17 are feasible choices for M. The minimum value of p* 
corresponding to M* = 15 can be computed as p* — 33 since the line M = 15 intersects L and C at 
(15, = 33) and (15, pi = 32.81), respectively. 



According to (29) any feasible s must satisfy the following inequality which was also derived in 

Wu HH] 

s 2 {t 2 - 2(f - to)n) - 2s{n -t)(t- t Q ) + (t- t ) 2 > 0. 
This in turn implies that 



s > 



{t-t Q )(n-t + y/n(n-d)) 
t 2 ~2n(t-t Q ) 



From (31 1 it also follows that a feasible value of s will exists only if 

i 2 - 2n(t - t ) > 0, 

which also leads to the same bound on the list decoding radius as derived in [TO] 

t < n — y n(n — d)) . 



Also from (31 1 we get the lower bound on s as 

{t-t )(n-t+ yjn{n-d)) 
Sl ~ i 2 - 2n(t - t ) 

Moreover, an upper bound on s was derived in |29j as 

t(2t Q - t) 
t 2 - 2n(t - t ) 



(30) 
(31) 

(32) 

(33) 

(34) 
(35) 



Thus any s, such that s; < s < s u will satisfy the condition (29). Now assume that for a particular 



s, the condition (29) is satisfied, i.e., L and C intersect at two different points (Mi,pi) and [M2.p1) 
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on the real plane. Without any loss of generality let us assume that Mi < M 2 . When L and C 
intersect at two different points on the real plane, there will exist a feasible solution if there is an 
integer M* such that M x < M* < M 2 , i.e., if 

LMiJ + 1< M 2 . (36) 

is satisfied, then any M* € [|_MJ + 1, \M 2 ~\ — 1] is a feasible choice of M. Now 



Clearly if (36 



according to (20), for a feasible choice of M = M* , it is desirable to find the minimum value of p 
so that U is minimized. Let the line M = M* intersect C and L at points (M* , pi) and (M*,ph) 
respectively. Since L intersects C from above, it must be the case that pi < ph- Although ph is a 
feasible choice for p, as used by Wu, we choose the minimum possible value as 

p* = [ Pl ] + 1. (37) 

We illustrate the method of computing the feasible values of the integer parameters, using a particular 
example, in Fig. [T] 

Now to find the optimal value of M and p such that MU is minimized, we need to compute p* 
and U* for all M* 6 [|_MiJ + 1, \M%\ — 1] and choose M* and p* that result in the minimum value 
of MU. We summarize the above procedure in Algorithm [6] that computes the values of minimum 
multiplicity s m ; n and the associated optimal z-degree M opt and weighted degree p op t- 



Algorithm 6 Compute optimal values of the integer parameters 
Input: n, k, t, fci, and 

Output: Minimum multiplicity s m i n and optimal z-degree M opt and weighted degree p pt- 
Compute w := k\ — k 2 , d := n — k + 1, to :— d/2. 

Initialize s :— max(s; = [_(< — t )(n — t + \Jn(n — d))/(t 2 — 2n(t — t ))\ + 1, 1) 

M opt := 00, p opt := 00, [/ opt = 00 
while no feasible solution is found do 
Compute N := ns(s + l)/2. 

if (ts - t + t Q ) 2 > A{N - ts)(t - t ) then 

(M 2 , Mx) := {(ts - t + 1 ) ± ^{ts - t + t ) 2 - 4(7V - ts)(t - t Q ))/2(t - t ) 
if [Mi\ + 1< M 2 then 

^min • S 

for M = [M X J + 1 to \M 2 ] - 1 do 

p := [N/{M + 1) + w/2M - 1J + 1 

U:= (p + l)(M + l) -w/2M(M + 1) 

if MU < M opt U opt then 

M opt := M, p opt := p, U opt = U 

end if 
end for 

return s min , M opt , p op t 
end if 
end if 
s := s + 1 
end while 



Complexity of Algorithm^ The complexity of the algorithm is dominated by the while loop and 
the for loop. Number of times the while loop is executed is bounded by s m i n . The for loop executes 

0(y/(ts - t + 1 ) 2 - 4(A - ts)(t - to)/(t - to)) = 0(ts) (38) 
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times. Moreover, the maximum list decoding radius is t = \n — \fn{n — d)) — 1] = 0(n). Thus 
Algorithm [6] computes the integer parameter values in time 0(ns 2 ). 

4.3 Computation of the message polynomial 

After constructing the bivariate polynomial, the solutions to the rational interpolation problem can 
be obtained by the rational factorization procedure of [55]. Clearly every solution (a, b) to the 
rational interpolation problem gives a valid error locator polynomial /' 2 ' = ag\^ + bg% ■ Given 
a valid error locator polynomial f^ 2 \ Wu's algorithm uses Forney's formula to compute the error 
magnitudes and hence the codeword. However, in our approach, the message polynomial can be 
computed in a simpler way: for every solution (a, b), it can be computed as 



m{x) = — 



ag^ + bgfy 



4. 4 Complexity 

We summarize the complete minimal list decoding algorithm in Algorithm [7] The computation of 

Algorithm 7 Minimal list decoding of (n, k) RS code using rational interpolation 
Input: Received word r = (r\, . . . , r n ). 

Output: A list of polynomials m of degree < k such that dn(c, r) is minimal, where c = 
(m(xi), . . .,m(x n )). 

1. Compute a minimal Grobner basis G = {<?i, 52} of M(r) with respect to the (0, k — l)-weighted 
top monomial order, with lpos((72) = 2 using Algorithm [3] (or using Algorithm [5] if re-encoding is 
used). Denote l\ := wdcgg! and £ 2 := wdegg 2 ; set j = 0. 

2. With t :— I2 — k + 1 + j, ki := £2 — h + j, and k 2 := j compute s min , M opt , and p opt using 
Algorithm [6j 

3. Construct Q(x, z) of wdeg 1 Q(x, z) = M opt and wdegj Q(x, z) = p opt passing through 
(xi, Zi)f =1 , with multiplicity s m i n using the Kotter algorithm from [21j . 

4. Compute all factors of Q(x, z) of the form z — a/b using the rational interpolation algorithm 
from [55] . 

5. If step 4 is successful, output all obtained quotient polynomials, i.e., polynomials m of the form 
m = 

-/ (1) // (2) ; Otherwise increase j by 1 and go to step 3. 

the minimal Grobner basis in step 1 using Algorithm [3] takes 0(n 2 ) operations. Algorithm [6] in 
step 2 takes 0{ns 2 ) time. The Kotter algorithm used in step 3 involves 0(MN 2 ) — 0(Mrrs 4 ) 
operations [H], where N is the number of constraints as defined in (l9| and M is the z-degree 
of the interpolating polynomial Q(x, z). The rational factorization in step 4 can be done in time 
0(n 3 / 2 s 7 / 2 ) [29]. Thus the overall complexity of the proposed algorithm is 0(MN 2 ). However, 
because of step 2, our list decoding algorithm optimizes MU. Since, more precisely, the Kotter 
algorithm involves memory complexity of O(MU) and time complexity of O(NMU), our algorithm 
uses less memory as well as computation as compared to Wu's method. The advantage of the 
proposed algorithm in terms of z-degree M and number of unknown coefficients U is illustrated in 
Example |4.1| 



Example 4.1 Consider the (127, 24) RS code defined over GF(2 7 ) with d = 104. Consider correct- 



ing t — 64 errors when ki = 15 and k 2 — 9. For this instance, Wu's algorithm using (22) computes 
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s = 2, which is also the minimum multiplicity. Now Wu's algorithm computes M = 5 and p = 72 
using (23) and (24), respectively. With these values, Wu's algorithm requires solving a system of 
N = 381 homogeneous equations involving U = 408 unknowns. In contrast, in our algorithm we 
find that when s min — 2, the line L intersect the curve C at points (3.3241, *) and (6.3426, *). Now 
for the feasible values of M* = 4, 5, 6, we get p* — 88, 78, 72 and U* — 385, 384, 385, respectively. 
Finally we get the optimal values as M opt = 4 and p op t = 88 with U op t = 385. 



5 Conclusions 

In this paper we have taken a parametric approach to the problem of minimal list decoding. The 
proposed algorithms have error correcting radius L, where L is the minimum of the Hamming 
distances between the received word and any codeword in C . There are several important features 
of the approach. Firstly, the minimality of L ensures that all solutions correspond to valid codewords 
and therefore we do not need to check for validity. The parameterization can also be used for general 
list decoding, however, then a check on the validity of the corresponding codewords needs to be 
carried out. Secondly, upon computation of a solution of the rational interpolation problem or, 
equivalently, of an error locator polynomial, we do not need to determine the error magnitudes via 
Forney's formula. Instead, solutions to the rational interpolation problem directly lead to message 
polynomials. Thirdly, we provide a geometric approach to optimize the integer parameters associated 
with the problem of rational interpolation. Since the interpolation step is the most computationally 
intensive task in list decoding, optimization of the integer parameters results in less computational 



as well as memory requirements. Finally, by using re-encoding as in sub-section 3.3 the approach 
lends itself well to the type of distributed source coding (DSC) proposed in [3]. 
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